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SOMINT:  AN  IMPROVED  MODEL  FOR 
STUDYING  CONDUCTING  OBJECTS  NEAR 
LOSSY  HALF-SPACES 


ABSTRACT 


This  report  summarizes  the  development  of  a moment-method  numerical  code  to  study 
wire  antennas  near  lossy  earth.  The  computation  time  of  this  code  is  greatly  reduced  by  two 
new  methods  for  evaluating  the  Sommerfeld  integrals.  In  the  evaluation  of  the  integrals  off 
the  earth,  a bivariate  interpolation  scheme  on  a prestored  grid  is  shown  to  give  very  quick 
and  accurate  results.  A new  series  solution  is  presented  to  evaluate  the  integral  for  values 
near  the  interface.  Test  results  of  this  program  are  given,  along  with  extensive  numerical 
data  on  vertical  half-rhombic  and  sloping  V antennas. 

f 

INTRODUCTION 


A problem  that  has  traditionally  resisted  a genuinely  practical  solution  despite  considerable  study  is  that 
of  modeling  conducting  structures  (antennas  or  scatterers)  located  near  a planar  interface  such  as  the  earth's 
surface.  Nearly  70  years  ago  Sommerfeld  (1909)  worked  out  the  basis  for  its  rigorous  solution  in  terms  of 
Fourier  integrals  of  cylindrical  wave  expansions.  The  Sommerfeld  integrals  have  been  extensively  studied,  and 
numerous  approximations  have  been  developed  for  them  for  various  combinations  of  the  problem's 
parameters.  1 However,  evaluation  of  the  Sommerfeld  integrals  for  the  ranges  of  source  and  observation 
points  that  typically  must  be  covered  for  any  self-consistent  description  of  the  structure's  current — an  integral 
equation  for  example — almost  inevitably  involves  some  parameter  combinations  to  which  such  approxima- 
tions do  not  apply.  Consequently,  the  only  feasible  approach  for  obtaining  a reliable  solution  in  such  cases  has 
been  to  numerically  integrate  the  Sommerfeld  integrals,  a computationally  inefficient  process. 

This  situation  has,  until  now,  kept  moment-method  models  of  wire  antennas  based  on  internal-equation 
formulations  too  expensive  for  routine  use  in  treating  the  half-space  problem.  In  this  report  we  present  a new 
approach  to  the  interface  problem  that  is  based  on  an  old  (interpolation)  numerical  procedure.  Our  interpola- 
tion approach,  which  we  refer  to  as  SOMINT  for  convenience,  is  shown  below  to  make  computer  modeling  of 
antennas  near  an  interface  now  affordable.  For  example,  relative  to  the  free-space  computation,  the  time 
penalty  for  rigorous  treatment  of  the  interface-reflected  fields  will  be  ultimately  a factor  of  only  5 to  10  vs  the 
former  value  of  100  to  1000.  Consequently,  a whole  new  range  of  problems  becomes  practicable,  including,  as 
illustrated  here,  radiation-pattern  control  through  ground-screen  design. 

The  second  section  of  this  report  provides  analytical  background  for  the  interface  problems  and  formula- 
tion of  the  integral  equations  upon  which  our  approach  is  based.  The  third  section  discusses  the  numerical 
procedures  employed,  including  the  interpolation  treatment  that  is  the  key  contribution  of  this  work. 
Numerical  results  are  given  in  the  fourth  section:  first  a series  of  test  cases  and  preliminary  calculation  and 
then  the  results  of  some  parametric  calculations  performed  for  two  Army  antennas,  the  vertical  half-rhombic 
and  sloping  V.  The  final  section  presents  our  conclusions  and  recommendations  for  future  work.  Two  appen- 
dices contain  representative  plots  of  radiation  patterns  and  current  magnitude  values  for  vertical  half-rhombic 
and  sloping  V antennas. 


METHODS  FOR  STUDYING  WIRE  ANTENNAS 

Preliminary  Considerations 

An  integral  equation  for  a wire  structure  can  be  derived  in  many  ways.  The  basic  procedure  involves 
writing  Maxwell’s  equations  in  integral  form  to  exoress  the  scattered  or  secondary  fields  in  terms  of  integrals 
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over  induced-source  distributions.  Expressing  the  secondary  field  over  loci  of  points  where  the  behavior  of  the 
total  field  (incident  or  primary  plus  secondary)  is  known  via  boundary  or  continuity  conditions  yields  an  in- 
tegral equation  for  the  induced  source  in  terms  of  the  primary  field.  Two  broad  general  classes  of  integral 
equations  result,  depending  on  whether  the  forcing  function  (primary  field)  is  electric  or  magnetic.  An  electric 
forcing  function  gives  rise  to  a Fredholm  integral  equation  of  the  first  kind  (so  called  because  the  unknown  ap- 
pears only  under  the  integral).  A magnetic  forcing  field  gives  rise  to  a Fredholm  integral  equation  of  the 
second  kind  (in  which  the  unknown  also  appears  outside  the  integral).  While  derivatives  of  the  unknown  may 
occur  as  well,  these  equations  are  commonly  termed  integral  equations,  rather  than  integro-differential  equa- 
tions. as  would  be  strictly  correct. 

Generally  speaking,  the  magnetic-field  type  of  integral  equation  has  been  found  better  suited  for  smooth, 
closed  surfaces  than  for  thin-plate  or  shell  geometries  and  wires.  2 The  converse  is  generally  true  for  the 
electric-field  type  of  integral  equation,  which,  accordingly,  is  most  commonly  employed  for  treating  wire 
structures.  Also  involved  in  developing  wire  integral  equations  are  three  approximations: 

• The  circumferential  current  is  negligible, 

• The  circumferential  variation  of  the  longitudinal  current  can  be  ignored;  and 

• The  thin-wire  or  reduced  kernel  can  be  used  in  place  of  the  actual  surface  integration. 

Many  analytically  equivalent  integral  equations  for  wires  based  upon  the  electric  field  can  be  derived. 
Three  of  the  most  commonly  employed  are  the  Hallen  or  vector  potential  type, 3 the  scalar-vector  (mixed) 
potential  version,  4 and  the  Pocklington  integral  equation.  5 All  are  solved  within  the  framework  of  the  mo- 
ment (or  matrix)  method,  but  each  exhibits  distinctive  characteristics  that  must  be  taken  into  account  in  its 
numerical  treatment.  The  Hallen  equation,  for  example,  can  produce  results  using  a pulse-current  basis  of  ac- 
curacy comparable  to  those  obtained  from  the  Pocklington  equation  solved  with  a three-term  (constant,  sine, 
and  cosine)  basis  for  simple  structures.  6 The  Hallen  equation  is  not,  however,  readily  extendable  to  the  com- 
plex geometries  that  the  Pocklington  equation  can  handle.  7 

Although  pulse-current 8'9  and  linear-current  10  bases  have  been  quite  widely  used,  and  can.  under 
suitable  circumstances,  be  essentially  equivalent,  they  are  not  as  efficient  for  modeling  traveling  wave  equa- 
tion, as  are  sinusoidal  bases  that  possess  nonconstant  derivatives  and  can  closely  resemble  the  actual  current 
solution.  Sinusoidal  bases  have  appeared  in  subsectional  or  subdomain  form  in  both  the  three-term  expansion 
mentioned  above  and  in  the  piecewise  sinusoidal 5 or  two-term  form.  Fourier  series  have  also  been  studied  as 
complete-domain  sinusoidal  bases,  but  have  not  been  widely  adopted  because  they  require  more  integration 
effort  than  subsectional  bases  and  can  lead  to  ill-conditioned  matrices.  8 

The  weight  or  test  functions  most  often  used  have  been  delta  functions,  although  Galerkin's  method  with 
both  linear  (two-term) 10  and  sinusoidal  (two-term) 5 functions  has  also  been  quite  widely  applied.  The  term 
“point  matching"  refers  to  the  use  of  delta-function  weights.  A comparison  of  numerical  convergence  rates  for 
several  common  methods  applied  to  a straight-wire  scatterer  is  shown  in  Fig.  1.  11 

In  addition  to  the  problem  of  choosing  basis  and  weight  functions,  other  special  aspects  of  the  numerical 
development  must  be  considered  when  selecting  a code  for  computer  modeling.  Three  of  these  aspects  are  dis- 
cussed below. 

Junction  Treatment 

Any  subsectional  approach  that  employs  either  finite-difference  operators  in  the  integral  equation  or  any 
other  basis  than  pulse-current  requires  special  consideration  of  both  simple  (two  wires)  and  multiple  (three  or 
more  wires)  junctions.  The  problem  essentially  is  to  relate,  in  some  physically  and  mathematically  reasonable 
way,  the  current  basis  of  each  subsection  (segment)  to  those  of  its  neighbors.  When  pulse  bases  are  used  in  the 
scalar-vector  potential  integral  equation,  the  finite-difference  operator  spans  two  segments  and  thus  leads  to  a 
charge  involving  the  two  corresponding  pulse-current  samples.  4 For  two-  or  three-term  bases,  the  condition 
of  current  amplitude  and  slope  continuity  at  each  simple  juction  leads  to  equations  that  permit  all  the  con- 
stants in  the  current  expansion  to  be  given  in  terms  of  current  samples  at  the  segment  junctions  or  centers.  6 
Yeh  and  Mei  12  developed  a slightly  different  handling  of  the  three-term  basis  in  which  the  current  is  ex- 
trapolated from  a given  segment  to  the  adjacent  segment  centers,  but  which  is  otherwise  basically  the  same. 

For  a multiple  junction,  the  treatment  can  become  considerably  more  involved.  The  pulse-basis  approach 
mentioned  above  has  been  extended  to  the  multiple  junction  9 by  dividing  the  total  junction  charge  between 
the  junction  segments  according  to  the  ratio  of  their  individual  areas  to  the  total  area  on  a logarithmic  basis  in 
radius  and  linear  in  segment  length.  This  approach  has  been  found  useful  for  the  three-term  expansion  as 
well.  6 The  two-term  expansions  have  been  applied  to  multiple  junctions  by  overlapping  M - 1 of  the  bases  a 
pair  at  a time  at  an  M-segment  junction.  10  Applications  of  the  three-term  expansion  to  the  multiple  junction 
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Fig.  I.  Convergence  rale  for  several  solution  methods  (Ref.  II). 


has  been  accomplished  by  MBAssociates  using  the  Yeh  and  Mei  simple-junction  procedure  by  incorporating  a 
composite  segment  having  the  averaged  length  and  total  current  of  the  M - 1 connected  segments.  13  A more 
elaborate  multiple-junction  approach  has  been  developed  for  the  three-term  expansion  by  Andreasen  and 
Harris.  14  Their  procedure  apparently  is  the  only  one  in  which  the  junction  geometry  plays  an  explicit  role  in 
establishing  the  current  relationships  at  the  junction.  Although  all  of  these  approaches  evidently  can  produce 
satisfactory  results,  there  is  little  or  no  direct  evidence  of  their  validity  in  terms  of  the  junction  current  and 
charge  distributions.  It  should  be  noted  that  the  numerical  results  have  been  found  in  some  cases  to  be  quite 
sensitive  to  the  junction  treatment.  6 Further,  the  above  list  by  no  means  exhausts  all  possible  alternatives  for 
the  junction  treatment. 

Source  Models 

Determination  of  quantities  such  as  absolute  gain,  efficiency,  radiated  power,  input  power,  etc.,  requires 
not  only  the  antenna  current  distribution  but  also  the  input  characteristics,  particularly  the  feedpoint  im- 
pedance (or  admittance).  The  feedpoint  admittance  can  be  found  in  various  ways,  but  using  the  integral- 
equation  approach  ordinarily  also  involves  defining  it  in  terms  of  source-region  current  per  unit  of  terminal 
voltage.  Calculation  of  this  quantity  requires  a realistic  source  model  that  not  only  provides  an  appropriate 
means  for  numerically  exciting  the  antenna,  but  also  permits  ready  evaluation  or  specification  of  the  effective 
terminal  voltage.  Thus  if,  as  in  a point-matching  procedure,  the  excitation  arises  as  a tangential  field  on  the 
source  segment  of  length  A,  the  driving  voltage  might  be  assumed  to  be  -E  'A  if  E 1 is  constant  on  the  source 
segment  and  zero  elsewhere.  This  assumption  may  be  invalid,  however,  with  the  result  that  the  actual  voltage 
can  be  obtained  only  by  integrating  the  tangential  field  in  the  vicinity  of  the  source  segment.  6 Somewhat  less 
ambiguity  should  arise  from  Galerkin-type  approaches  where  the  boundary  conditions  are  integrated,  so  that 
the  classical  delta-function  source  might  be  numerically  approximated.  An  alternative  source  model  for  point- 
matching is  provided  as  a current-slope  discontinuity,  which  also  approximates  a delta-function  source  field. 
The  current  bases,  junction  treatment,  and  weight  functions  can  all  influence  the  usefulness  of  these  alter- 
native source  models.  In  case  of  uncertainty,  once  the  current  distribution  has  been  found,  the  impedance  can 
be  computed  from  the  classical  emf  method,  although  at  the  expense  of  the  additional  integration  this  entails. 
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Integration 

Integration  is  understandably  an  essential  part  of  the  moment  method;  it  is  involved  in  applying  the  in- 
tegral operator  to  the  current  bases  and,  in  a Galerkin  method,  evaluating  the  inner  product  of  these  results 
with  the  weight  functions.  For  most  wire  programs,  these  operations,  which  lead  to  the  generalized  impedance 
matrix,  dominate  the  total  solution  time  for  numbers  of  unknowns  less  than  ~200.  It  is  thus  important  that 
the  integration  time  be  minimized  consistent  with  the  overall  accuracy  requirements. 

One  way  to  approach  this  goal  is  to  choose  appropriate  bases  and  weight  functions.  The  two-term 
sinusoidal  current  basis,  for  example,  requires  no  numerical  integration  when  the  Pocklington  integral  equa- 
tion is  used  together  with  point  matching.  This  particular  combination  is  not  very  accurate,  however.  11 
Adding  the  constant  term  leads  to  much  better  results,  with  the  slight  additional  expense  of  the  numerical  in- 
tegration required  to  find  the  longitudinal  field  of  this  current  term;  the  radial  component  can  be  analytically 
expressed.  Alternatively,  use  of  a sinusoidal  weight  function  5 also  gives  much  improved  results  and  sur- 
prisingly requires  numerical  integration,  at  most,  of  sine  and  cosine  integrals.  The  piecewise  linear  basis  used 
with  the  scalar-vector  potential  equation  cannot  be  analytically  integrated,  but  good  results  are  obtained  with 
four-point  rectangular  integration  of  both  the  operator  and  inner-product  integrals.  In  addition,  instead  of  ap- 
plying a numerical  integration  of  the  self-term,  a series  expansion  that  gives  a closed-form  expression  is 
used.  15  When  numerical  integration  is  used,  various  adaptive  routines  and  special  techniques  are  available  to 
improve  efficiency.  16 

Wire  Antenna  Analysis 

Examining  the  many  options  available  to  the  analyst  concerning  the  integral  equation  to  be  selected  and 
its  numerical  treatment  in  developing  a computer  model  for  application  to  wire  antennas  would  take  this  dis- 
cussion beyond  reasonable  lengths.  Accordingly  our  attention  is  primarily  directed  to  an  approach  based  on 
the  Pocklington  integral  equation  solved  using  a three-term  subsectional  basis  (constant,  sine,  and  cosine)  and 
point-matching.  Unless  otherwise  indicated,  antenna  sources  are  introduced  as  tangential  electric  fields,  with 
the  Yeh  and  Mei 12  form  of  current  extrapolation  used  for  simple  junctions  and  the  Curtis  method  for  multiple 
junctions.  Both  the  source  model  and  junction  (simple  and  multiple)  treatment  used  in  this  code  may  exhibit 
deficiencies,  but  when  applied  with  care  (e.g.,  equal  segment  lengths  near  sources  and  at  multiple  junctions) 
the  code  has  proven  valid  and  reliable.  A brief  overview  of  the  relevant  equations  and  numerical  treatment 
used  for  free  space  and  various  interface  theories  and  some  special  topics  is  given  in  this  section.  Numerical 
results  appear  in  a later  section. 

Infinite,  Homogeneous,  Isotropic  Media 

The  Pocklington-type  integral  equation  for  a wire  structure  of  contour  C(r)  can  be  expressed  in  the  form 
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where,  as  usual,  k is  the  inifinite-medium  wave  member,  the  permeability  and  permittivity  are  denoted  by  n 
and  (.  a(r)  is  the  wire  radius  at?,  and  E 1 is  the  incident  electric  field. 

Reduction  of  this  equation  to  matrix  form  involves  these  seven  steps: 

(1)  Approximating  C(r)  as  a piecewise  linear  sequence  of  N segments  of  length  Aj  (i  - 1 N),  so  that 

N 

C(r)  - £ A&  , 
i=l 

with  Sj  the  unit  tangent  vector  to  C(r)  at  r = (use  of  straight  segments  is  not  mandatory,  but  it  is  very  con- 
venient in  simplifying  the  current  integration); 

(2)  Introducing  the  subsectional  bases 

l,(s')  = Aj  + Bj  sin[k(s'  - Sj)]  + Cj  cos[k(s’  - Sj)] 

to  represent  the  unknown  current  (the  final  unknowns  will  be  the  N sampled  current  values  1;  = A;  + Cj  [t  = I , 
N)  at  the  center  of  each  of  the  N segments); 

(3)  A current  interpolation  procedure  whereby  the  individual  Ar  Bj,  and  Cs  constants  are  expressed  in 
terms  of  the  sampled  current  values; 

(4)  Use  of  the  N delta-function  weights  S(s  — Sj)  (j  = 1 N)  to  obtain  an  N -order  impedance  matrix 

of  N independent  field  equations  (note  that  the  weight  functions  sample  the  field  at  the  segment  centers,  and 

are  thus  “collocated"  with  the  current  sample  locations);  A . 

(5)  Specification  of  the  N incident  or  primary  field  vector  components  Ej  - E (sj)  • Sj  (j  - 1,  N), 

which  are  the  tangential  fields  at  the  N segment  centers; 

(6)  Matrix  manipulation  to  obtain  an  admittance  equivalent  of  the  impendance  matrix;  and 

(7)  Computation  of  the  current  distribution  and  whatever  field  components,  if  any,  are  desired. 

The  total  computer  solution  time  is  well  approximated  by  AN  2 + BN  \ where  A corresponds  to  step  (4) 

and  B to  step  (6).  For  the  code  under  consideration  here  and  for  a CDC  7600  computer,  A«4X  10  s and  B 
« 2 X I0-6  s. 

Perfectly  Conducting  Half-Space 

As  written,  Eq.  (1)  applies  to  wire  structures  excited  as  antennas  or  scattered  and  located  in  infinite, 
isotropic,  homogeneous  media  of  arbitrary  (possibly  lossy)  permittivity  and  permeability.  It  can  easily  be  ex- 
tended to  permit  the  modeling  of  magnetic  or  electric  image  planes.  For  example,^the  perfectly  conducting 
ground  analog  of  Eq.  (1)  is,  for  an  antenna  elevated  above  a ground  plane  at  z = 0, 

§ • E!(s)  = ^ l(s')  [G0(sy)  + G|(s,s'  )]ds'  , W 

C(r) 


where 


r*  = if  - r*i  , 


5 


r 


T*(x,y,z)  = 7'(x,y,-z)  , 

A.*  = VC(f'*) 

~ lVC(r'*)| 

Similar  forms  can  be  written  for  a magnetic  interface  and  for  an  interior  right-angle  corner.  If  the  corner 
angle  is  otherwise  arbitrary  but  related  to  n as  an  integer  multiple,  a discrete  spectrum  of  angular  images  is  ob- 
tained, but  the  essence  of  the  integral  equation  form  is  preserved.  Precisely  the  same  line  of  approach  can  also 
be  used  for  interior  problems  where  the  wire  structure  is  located  between  two  parallel  magnetic  or  electric 
planes.  17 

Imperfectly  Conducting  Half-Space 

A prob'em  that  is  not  so  computationally  simple,  but  that  is  of  perhaps  greater  practical  interest,  is  that 
of  an  antenna  located  (buried  or  elevated)  near  the  ground-air  interface.  This  is  a topic  of  considerable 
longevity  in  electromagnetics;  a formal  solution  to  the  problem  was  given  in  1909  by  Sommerfeld.  18  The 
numerical  complexity  of  evaluating  the  Sommerfeld  integrals  (which  appear  in  the  integral  equation  kernel) 
for  arbitrary  source  and  observation-point  locations  and  ground  parameters,  however,  has  prevented  the 
Sommerfeld  theory  from  being  routinely  used  for  such  problems.  Consequently,  while  some  progress  has  been 
made  in  applying  the  Sommerfeld  theory,  alternative  approaches  to  the  antenna-ground  problem  have  also 
been  pursued.  These  various  methods  are  briefly  discussed  below. 

The  Sommerfeld  Theory 

Details  of  the  steps  in  deriving  the  Sommerfeld  integrals  may  be  found  elsewhere.  19  Here  we  simply  write 
one  version  of  Eq.  (1)  that  accounts  for  the  interface  reflected  field  via  the  Sommerfeld  theory;  alternative 
forms  are  also  available,  differing  essentially  in  how  the  perfert-ground  image  terms  are  handled.  It  is 


/ l(s’)ds'  X <GQ  (s,s')  + G,  (s,s'*> 
■'C(T)  l 


+ ^2  ak)  sin(5'  ®h  7.  ~ C0V  gVz 


sin0'  jsin|3  cosfa  - «')  + ± j^]  gN  ,1 , 


(3) 


where  a = a(r)  and  0 = 0(r)  are  the  direction  angles  of  the  wire  at  r,  t'  is  the  horizontal  projection  of  s',  Jn  is 
the  Bessel  function  of  order  n,  and 


gHt  2 


r~L 

l yo  + 


y . 0 


— y n(z+z') 

J„(Ap)e  0 dX  , 


g„  = -cos(0  - «')  f -2S. Ii_  x fj  (*p)e  7°<Z+?-  ’ X2dx]  , 

4,  koT'l  + k^0  L 
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g Vz 


r 


+ k 


I TO 


p = V(x  - x')2  + (y  - y')2  + a2  , 
0 = tan- 1 [(y  - y')/(x  - x')]  , 


T0  = V x2  " ko  > 

T,  = V x2  - kf  ■ 

where  ko  and  k)  are  the  propagation  constant  above  and  below  the  interface,  respectively. 

The  presence  of  the  double  integral  in  Eq.  (3),  particularly  the  Sommerfeld  portion,  makes  it  quite  time- 
consuming  and  sensitive  to  evaluate.  In  spite  of  that,  the  basic  moment  method  can  be  used  to  solve  it,  but,  in 
addition  to  the  usual  constraints  imposed  on  current  sampling,  it  is  necessary  to  take  into  account  the  source 
distance  from  the  interface. 


Series  Solution  to  Sommerfeid  Integrals 

In  the  study  of  antennas  over  lossy  interfaces,  the  Sommerfeld  integrals  need  a complex  integral  that  has 
an  infinite  path  of  integration.  There  are  three  methods  to  evaluate  these  integrals.  The  first  procedure  is  the 
asymptotic  expansion  of  the  integral  developed  by  Norton.  20  In  this  expansion  the  fields  are  written  in  terms 
of  the  inverse  powers  of  the  field-point-source-point  distance.  The  asymptotic  expansions  are  accurate  only 
for  large  field-point-source-point  distances.  The  quasistatic  method  presented  by  Banos  1 is  an  approximation 
that  sets  k|  = 0 in  the  integrand;  the  modification  integrals  are  then  evaluated  in  closed  form.  This  procedure 
is  supposed  to  work  when  the  field  points  and  source  points  are  very  close,  but  the  accuracy  of  this  method  has 
I never  been  established.  The  last  method  is  the  numerical  integration  of  the  Sommerfeld  integrals.  A numerical 

code  using  this  procedure  has  been  developed  by  Lytle  and  Lager. 21  This  procedure  gives  accurate  results  for 
i most  values  of  field-point-source-point  distance,  but  is  is  a very  slow  numerical  procedure. 

A general  numerical  procedure  to  evaluate  the  Sommerfeld  integral  for  all  values  of  field-point-source- 
point  distance  is  discussed  in  a later  section  of  this  report.  In  this  section  we  present  a new  series  solution  that 
t is  later  used  to  evaluate  the  Sommerfeld  integral  when  the  source  point  and  field  point  are  near  the  interface. 

P The  fields  from  a source  point  and  field  point  on  the  free-space  side  of  the  interface  can  be  written  in  terms  of 

i two  integrals: 


— TT-  e" V Jo<0x> xdx  - <4a> 

70  T| 


r“'  '.**> » «■» 
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An  oull'ne  ot  the  procedure  to  find  the  series  expansion  for  the  former  integral  is  presented  below. 
The  development  of  a new  series  representation  begins  by  writing  y0  = (X  2 - k|  + >1')  >/2  in 

1 

>0  + >1 


and 


-7 

e 


0Z 


where  'L  = kf  - kg.  Then  these  new  terms  are  expanded  by  Taylor  series  in  The  first  Taylor  series  gives 


1 

70  + T, 


2 G„,  [k;  - k«]m 


(5) 


when  the  first  coefficient  Gq  - 1/2.  The  remaining  G’s  can  be  found  by  a recursion  procedure.  Let  the  n Ih 
term  be  expressible  as 


c„  ■ £ C 

m=l 


(6) 


then  the  n+  1 lh  term  is 
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where 
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The  first  term  in  the  above  recursion  is 


g!  8 


The  Taylor  series  for  the  exponential  term  gives 

i c.  (ki  - kiT  • 


-r0z 

e u = e 


m=0 


(9) 


where Gq  = 1 . The  remaining  G’s  can  be  found  by  a recursion  procedure.  Let  the  n th  term  be  expressed  as 


p=i  i 


2m-p 


GO) 


then  the  n + 1 lh  term  is 


n+  I 


n+1 

= 2 

p=i 


gn+1 
P _.2n- 


■in-p 
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where 


*7H 


.2"  -.L  gn 

2(n  + 1) 


(12a) 


Dn+  I = _ J p"  _ 2t1  ~ P o" 

*P  2(n  + 1)  *P-l  2(n  + I)  *P 


fP  = 2,  3 n 

\n  = 2,  3,  . . n 


(12b) 
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The  first  term  in  the  recursion  is 


The  Taylor  series  represented  in  Eqs.  (5)  and  (9)  are  then  multiplied  together  to  find  a single  expansion  for  the 
term  e _‘>'0  */  + To  *n  the  integrand  of  Eq.  (4a).  This  new  representation  of  the  integrand  is  an  infinite  series. 

When  the  series  and  integration  are  interchanged,  the  integral  in  Eq.  (4a)  can  be  expressed  as 
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p=l 

The  term  Lm  in  Eq.  ( 14)  is  the  integral 


2 / J..  (Xp)  XdX  . 

I 7 m o' 
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To  evaluate  the  series  in  Eq.  (14),  we  must  first  find  a method  to  evaluate  the  above  integrals.  To  accomplish 
this  we  now  introduce  an  operator 


r = — — 

k , Sk,  ' 


When  this  operator  is  used  on  L,  and  L2  it  can  be  shown  that 


T T-1 

. j [r!  ” L,]  - ]T  KJ'z2"-'  [r2'r-  L,j  , 


...  ■-  i Hr 


f’  ILL  - 
r2  "P  L 


,J  + ^ KTz2,p-,)  f2 +' P 


where  m is  an  even  integer.  In  the  above  expression,  V m means  that  the  operator  T is  applied  m times  to  L:  or 
L|. 

The  integral  L|  is 


ik,(p2+z2,'/2 


L,  = 2 -- 


(p2  + z2)l/2 


By  applying  the+‘ m operator  to  L(.  we  can  show  that 
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The  I s in  the  above  equation  can  be  generated  from  a recursion  relationship. 
The  first  two  I's  are 

I,1  = (-i) 

and 

\\  = (02[(k,P)2  + (k,z)2]1/2  . 

The  remaining  term  can  be  generated  for 

in+l  _ in 
‘l  ~ *1  ’ 

in+I  _ in 
*2  l2  ’ 


(19a) 


(19b) 
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for  n > 3 and  p = 3,  4,  ....  n , 
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Now  that  we  have  a method  to  evaluate  the  f nL),  we  can  turn  our  attention  to  the  other  terms  that  still  need 
evaluation. 

The  integral  L2  can  be  shown  to  be 


? [Z  e^1  (P2  + z-2)^2 
‘ l (P2  ♦ Z'2)1/2 


dz'  + itrH},1*  (kjP)  . 


(20) 


When  the  F operator  is  applied  to  the  above  expression  one  finds 
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The  term  H,,*1*  is  the  Hankel  function  of  first  kind  in  the  order.  To  find  the  integer  in  Eq.  (21)  we  expand  the 
integrand  in  terms  of  z'  and  then  perform  the  integrand  analytic.  Then  •Lm  operates  on  the  new  series 
representation  for  the  integral. 


NUMERICAL  METHODS 


In  this  section  we  discuss  the  numerical  method  used  to  evaluate  Sommerfeld  integrals  in  this  study.  This 
approach  is  unique  in  that,  when  the  field  point  and  source  point  are  above  an  interface,  all  of  the  electrical 
field  from  the  Sommerfeld  integrals  can  be  written  in  terms  of  two-dimensional  space.  To  prove  this,  we  need 
only  look  at  the  Sommerfeld  integral 
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To  + T'l 


-7n(z+z’) 

e 0 Jp(Xp)XdX 


(22) 


where  p = \(x  - x') 2 + (y  - y')2.  The  points  (x,y,z)  and  (x\  y’,  z')  are  the  field  point  and  source  points, 

respectively.  All  values  of  the  above  integral  can  be  described  in  terms  of  two  variables  p and  z + z'.  We  now 
introduce  a new  two-dimensional  space  in  p and  z + z';  it  is  the  solution  space  for  Sommerfeld  integrals.  The 
solution  space  is  shown  in  Fig.  2. 

In  the  framework  of  the  solution  space  we  must  now  evaluate  the  Sommerfeld  integrals.  When  the  point 
in  the  solution  space  is  described  by  1.0  < p/\0  or  1.0  < z + z'/X0  the  asymptotic  series  developed  by 
Norton  20  is  used  to  evaluate  the  fields.  The  asymptotic  series  is  an  expansion  in  inverse  powers  of  the  field- 
point-source-point  distance:  therefore,  it  provides  a good  method  to  evaluate  the  integrals  when  the  distances 
are  large.  The  asymptotic  methods  do  not  work  well  when  the  field  point  and  source  point  are  less  than  one 
wavelength.  Figure  2 depicts  the  region  in  the  solution  space  where  the  asymptotic  method  is  used  to  evaluate 
the  Sommerfeld  integrals. 

To  gain  some  insight  into  the  method  used  for  evaluating  the  Sommerfeld  integrals  in  other  regions  of  the 
solution  space,  we  evaluated  the  various  components  of  the  electric  field  for  the  horizontal  and  vertical  dipoles 
in  the  solution  space.  The  evaluation  of  the  Sommerfeld  integrals  was  obtained  by  using  a code  developed  by 
Lytle  and  Lager 21  that  performs  the  integration  numerically.  A three-dimensional  plot  in  the  solution  space  of 
the  real-part  of  the  electrical  field  produced  by  a vertical  dipole  is  shown  in  Fig.  3,  which  plots  the  value  of  the 
function  in  terms  of  p,  z + z’.  The  observation  that  the  electric  field  value  plotted  in  Fig.  3 is  rather  well- 
behaved,  except  near  the  origin,  led  us  to  a procedure  known  as  the  “bivariate  interpolation  of  the  electrical 
fields.”  This  procedure  begins  with  division  of  the  solution  space  into  a set  of  rectangular  grids.  The  electrical 
field  for  a given  f,  tr,  and  a is  found  at  each  grid  point  in  the  solution  space  prior  to  evaluating  a given  antenna. 
W'hen  the  moment-method  code  needs  a field  value  not  on  the  grid,  a four-point  bivariate  interpolation  is  used 
to  evaluate  the  field  from  the  four  known  field  values  at  the  corners  of  the  grid.  This  procedure  is  outlined  in 
Fig.  4.  To  demonstrate  the  accuracy  of  this  approach,  we  show  in  Fig.  5 a three-dimensional  plot  of  the  error 
function  when  the  bivariate  interpolation  is  used  to  find  the  electrical  field  plotted  in  Fig.  3.  In  this  test,  a grid- 
spacing  of  0.5X  produced  a 10%  agreement.  In  general,  smaller  grid  spacing  produces  better  agreement. 

The  bivariate  interpolation  of  the  electrical  fields  provides  good  answers  in  the  solution  space  for  values 
ofp  and  z + z'  in  the  two  regions  where  [0. 1 < p/X0  < 1 .0  and  0 < z + z'/X0  < I ] and  [0  < p/A0  < 0. 1 and  0. 1 
< z + z'/Xq  < 1.0).  The  size  of  the  grids  used  in  this  calculation  is  shown  in  Fig.  6.  Bivariate  interpolation  re- 
quires a computation  time  of  about  100  ps  on  a CDC  7600  for  a single  set  of  field  values,  while  the  code  that 
numerically  integrates  the  Sommerfeld  integrals  correspondingly  requires  0.1  s.  (Note  that  after  the  solution- 
space  grid  has  been  evaluated  for  given  f,  «r,  and  a,  these  values  can  be  saved  for  future  use;  they  are  indepen- 
dent of  the  type  of  antenna  studied.)  The  electric  field  determined  from  the  Sommerfeld  integrals  has 
singularities  when  both  p = 0 and  z + z'  = 0.  The  series  representation  of  the  Sommerfeld  integrals  discussed 
earlier  isolates  these  singularities.  Therefore,  the  series  representation  is  used  to  evaluate  the  field  near  the 
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Fig.  2.  Two-dimensional  solution  space  used  to  evaluate  the  Som- 
merfeld  integrals.  The  various  methods  of  evaluating  the  integrals 
and  their  usable  regions  are  also  shown. 


Fig.  4.  The  four-point,  bivariate  interpolation  procedure.  The 
function  f is  found  in  terms  of  its  four  known  end-point  values  f0  0, 

f l,*  f«.|.  »nd  fM.  f = (1  - pH  I - q>f0.o  + P<<  - + 1<i  - 
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Fig.  3.  Three-dimensional  aspect  plot  of  Sommerfeld-integral 
portion  of  radial  electric  field  for  vertical  dipole.  The  real  part  of 
the  electric  field  is  plotted  on  the  vertical  axis  as  a function  of  z + 
z'  and  p.  The  electric  parameters  are  f = 10  MHz,  i,  = 9.0,  a = 
10  2 mho/m,  0.1  A < p < 2.0  X,  and  0.1  X < (z  + z')  < 2.0  X.  The 
maximum  and  minimum  field  values  are  also  shown. 

I error  lmax  = 0. 102 


yT7 V 


X 

Fig.  S.  Aspect  plot  of  difference  between  the  Sommerfeld-integral 
field  values  and  interpolated  values  for  real  part  of  radial  electric 
field  of  vertical  electric  dipole.  The  electric  parameters  are  the 
same  as  those  in  Fig.  3.  The  separation  between  interpolation  grid 
points  is  0.5  X.  The  interpolation  was  performed  for  the  region 
defined  by  0.1  X < (z  + z')  < 2.0  X,  excluding  the  small  region  0.1 
X < p < 0.5  X and  0.1  X < (i  T i’|  « 0.5  X. 
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Kig.  6.  The  grid  sizes  used  in  the  solution  space.  The  appropriate 
grid  sizes  are  shown  in  their  usable  regions  as  indicated. 


Table  1.  Comparison  of  the  magnitude  of  electric 
fields  evaluated  by  the  series  solution  with  those 
evaluated  by  numerical  integration  of  the  integral; 
f = MHz,  er  = 9,  and  a = 10-3  mho/m.  The  two 
fields  shown  are  the  radial  and  z-components  from 
a vertical  dipole. 
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origin  in  the  solution  space.  In  the  numerical  code  used  to  evaluate  this  field,  the  first  ten  series  terms  are  used. 
T able  I shows  a comparison  of  electrical  field  magnitudes  evaluated  by  the  series  method  with  those  generated 
by  numerical  integration.  From  Table  1,  it  is  obvious  that  the  series  solution  works  only  when  0 < p/X0  < 0.01 
and  0 < z + z'/X0  < 0.01.  The  CDC  7600  computer  time  required  to  evaluate  the  series  is  about  0.01  s while 
the  numerical  integration  requires  0. 1 s.  The  region  in  the  solution  space  occupied  by  the  series  solution  is 
shown  in  Fig.  2. 

As  Fig.  2 shows,  a region  exists  where  no  reliable  method  is  available  to  evaluate  the  Sommerfeld  integral. 
It  lies  between  the  regions  where  the  bivariate  interpolation  on  the  electrical  field  works  and  where  the  series 
solution  works.  In  this  region,  some  of  the  electric  field  may  vary  by  two  orders  of  magnitude.  Therefore,  a 
straight  bivariate  interpolation  on  the  electrical  field  is  inappropriate.  To  overcome  this  difficulty  we  in- 
troduce the  quantity 

R = ^ • (23) 

se 

where  the  Eex  and  Ese  are  the  electrical  fields  from  the  numerical  integration  and  series,  respectively.  The  new 
parameter  R is  a slowly  varying  function  in  the  region  under  consideration.  Since  the  quantity  R is  a well- 
behaved  function,  we  can  use  the  four-point  bivariate  interpolation  on  R.  This  procedure  consists  of  storing 
the  function  R on  a grid  of  points  in  the  solution  space.  This  grid  of  points  is  filled  before  any  moment-method 
code  is  used  to  study  the  antenna.  When  the  moment-method  code  needs  values  of  the  electric  field  at  a point 
not  on  the  grid  points,  a four-point  bivariate  interpolation  is  used  to  find  the  function  R at  this  point.  Thenjhe 
series  representation  of  E^,  calculated  at  the  point  off  the  grid,  is  multiplied  by  the  interpolated  value  of  R to 
obtain  an  interpolated  value  for  Eex.  We  call  the  procedure  “bivariate  interpolation  on  the  parameter  R.“  The 
CDC  7600  computer  time  required  to  evaluate  method  is  about  0.0 1 s.  The  grid  size  used  in  this  calculation  is 
shown  in  Fig.  6.  (Note  that,  for  a fixed  f,  «r,  and  <r,  the  value  of  R is  independent  of  the  antenna's  geometry. 
Once  the  values  of  R on  the  grid  point  for  a fixed  f,  tr,  and  a have  been  evaluated,  they  can  be  stored  on 
magnetic  tape  for  future  use.  This  saves  the  grid  fill  time  in  the  study  of  other  antennas.) 

TEST  CASES  AND  RESULTS 

Three  types  of  calculations  are  summarized  in  this  section.  The  first  type  involves  test  cases  to  a:sess 
operation  of  the  overall  program;  it  is  intended  to  provide  a check  on  numerical  accuracy.  The  second,  by  con- 
trast, is  intended  to  obtain  some  preliminary  results  prior  to  performing  parametric  calculations  on  the  an- 
tennas of  interest,  which  are  covered  in  the  third  section. 
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Test  Cases 


One  of  the  greatest  uncertainties  associated  with  using  a new  computer  model  is  validating  its  results. 
Thus,  much  of  our  attention  during  the  early  stages  of  this  work  was  focused  on  test  calculations  to  check  the 
Sommerfeld-integral  routine,  the  interpolation  procedure,  and  the  overall  computation  itself.  Calculations 
pertaining  to  the  first  two  cases  are  reported  in  the  previous  section.  Here  we  discuss  results  obtained  in  con- 
nection with  the  overall  computational  accuracy. 

Ideally,  experimental  results  provide  probably  the  best  and  most  convincing  data  for  use  in  validating  a 
computer  model.  Realistically,  however,  such  data  are  neither  extensive  nor  readily  available;  those  that  are 
obtainable  are  of  limited  scope  and  uncertain  applicability.  Therefore,  while  our  ultimate  goal  remains  that  of 
validating  the  computer  model  using  experimental  data,  the  test  cases  presented  here  are,  instead,  of  a 
numerical  nature. 

Several  kinds  of  numerical  tests  might  be  performed  to  check  the  computer  model's  accuracy.  One  of  the 
more  obvious  is  simply  to  examine  the  trend  of  the  numerical  results  for  obviously  incorrect  behavior  as  a 
function  of  the  various  parameters.  An  example  of  this  is  shown  in  the  preceding  section,  where  the  inter- 
polated fields  are  found  to  exhibit  anomalous  behavior  along  the  z + z'  = 0 and  />  = 0 lines.  This  kind  of  test 
may  be  regarded  generally  as  a necessary,  but  not  sufficient,  condition  to  be  met  by  a valid  calculation.  It  can 
be  especially  useful  in  the  preliminary  stages  of  code  development  or  as  a diagnostic  w hen  judgmentally  in- 
valid results  are  obtained  as  the  overall  computational  output. 

Another  useful  numerical  test  is  to  compare  results  from  the  new  code  with  those  obtained  from  the 
earlier  Lawrence  Livermore  Laboratory  (LI.L)  Sommerfeld  code,  which  did  not  use  interpolation.  This  was 
done  for  several  cases,  for  example,  horizontal  dipole  and  Beverage  antennas.  SOMINT  produced  results  that 
generally  are  numerically  indistinguishable  from  these  alternate  calculations  at  LLL,  21  as  shown  in  Figs.  7-10. 

Numerical  tests  can  also  be  based  on  limiting  cases  of  the  associated  parameters  where  the  answers  are 
known  analytically.  For  example,  when  the  relative  permittivity  and  conductivity  of  the  ground  approach 
their  free-space  values,  the  antenna  current  should  be  independent  of  height  above  the  interface.  The 
Sommerfeld-integral  fields  must  then  also  exactly  cancel  the  perfect-ground  image  fields,  so  that  only  the  free- 
space  contributions  remain.  It  is  convenient  to  check  this  by  examining  the  input  impedance  of  an  antenna  as 
it  is  lowered  towards  the  now  fictitious  interface.  Results  obtained  from  this  kind  of  calculation  are  shown  in 
Fig.  II.  Such  tests  proved  very  useful  in  correcting  numerical  inconsistencies  originally  present  in  the 
program,  especially  those  that  were  important  only  where  z + z'  -*  0. 


Preliminary  Results 

Antennas  operated  in  the  HF  range  and  located  near  the  earth's  surface  generally  are  excited  and/or  ter- 
minated with  respect  to  ground  potential.  This  may  be  accomplished  in  several  ways;  for  example,  a ground 
stake,  counterpoise,  ground  screen,  or  ground-return  wire  may  be  used.  The  first  two  methods  involve  wires 
that  penetrate  the  interface,  the  treatment  of  which  is  somewhat  uncertain.  In  addition,  modeling  them  re- 
quires computing  the  fields  transmitted  across  the  interface,  to  which  the  interpolation  procedure  developed 
here  does  not  directly  apply. 

Our  approach  has  thus  been  to  concentrate  performing  computations  for  antenna  models  whose  wires  lie 
entirely  on  the  free-space  side  of  the  interface.  Given  this  constraint,  we  can  model  configurations  which  in- 
corporate ground-return  wires  or  ground  screens  so  long  as  their  height  h above  the  ground  equals  or  exceeds 
that  previously  shown  to  be  necessary  (Fig.  1 1).  As  observed  in  connection  with  Fig.  1 1,  the  minimum  height 
for  valid  results  (<3  X 10  _4X0)  is  small  enough  that,  for  practical  purposes,  the  wires  are  at  the  interface  * 

One  other  situation  can  also  be  modeled  (albeit  only  approximately):  the  ground  stake  excitation /ter- 
mination mode.  In  this  case,  the  antenna  wire  is  terminated  at  the  interface  and  the  effect  of  the  ground  stake 
is  approximated  in  two  ways.  First,  the  current  at  the  end  of  the  wire  is  constrained  to  have  zero  slope,  a con- 
dition rigorously  correct  only  for  a perfect  ground,  but  a behavior  approximately  simulated  by  an  actual 
ground  stake.  Second,  the  fields  due  to  the  ground-stake  current  are  ignored. 

* Besides  the  limit  imposed  by  Sommerfeld-field  accuracy,  there  is  also  a minimum-height  limit  due  to  the  thin-wire  approximation.  This 
implies  that  for  h<,IOa.  the  problem  must  be  reformulated  if  questionable  results  are  to  be  avoided.  Of  course,  if  a hare  wire  is  lowered  to 
h = a and  thus  comes  into  contact  with  a possibly  lossy  ground,  still  further  reformulation  is  required  to  account  for  their  ohmic  contact 
These  aspects  of  the  problem  are  beyond  the  scope  of  the  present  work. 
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Fig.  7.  Three  antennas  used  to  compare  the  numerical  performance  of  the  interpolation  approach  with  the  orginial  Sommerfeld  mode:  (a) 
slanted  dipole  above  interface;  (b)  horizontal  wire  0.01  A above  interface  loaded  near  the  end;  (c)  horizontal  wire  0.05  A above  interface  loaded 
near  the  end.  The  current  distributions  are  presented  in  the  following  figures.  Small  differences  are  observable  between  results  obtained  from 
the  two  approaches,  but  the  basic  behavior  is  the  same. 

We  performed  calculations  pertaining  to  these  issues  for  monopole  and  vertical  half-rhombic  antennas. 
Some  of  the  results  are  presented  in  Figs.  12-22.  Based  on  the  outcome  of  these  preliminary  calculations,  we 
chose  the  ground-stake  approximation  to  model  the  grounded  ends  of  the  antennas  selected  for  parametric 
study.  At  this  point  of  development  in  the  program,  it  represents  the  least  certain  element  in  the  calculations. 


Parametric  Results 

We  performed  parametric  calculations  to  evaluate  ground-screen  effects  and  to  characterize  two  antenna 
types,  the  slanted  V and  vertical  half-rhombic.  The  results  obtained  are  summarized  below. 

Ground-Screen  Effects 

Ground  screens  can  provide  a terminal  against  which  to  drive  or  terminate  an  antenna.  More  impor- 
tantly, they  can  provide  impedance  and  pattern  control.  Impedance  is  controlled  hen  variations  due  to 
changes  in  local  ground  conditions  can  be  reduced  or  eliminated,  permitting  the  feedpoint  current  to  become  a 
function  of  the  screen  parameters.  Pattern  control  is  possible  because  the  ground-screen  currents  add  their 
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Position  — wavelengths 


t ig.  8.  Magnitude  of  (he  current  on  the  antenna  shown  in  Kig.  7a. 
I he  interpolation  solution  to  the  problem  is  compared  to  the  exact 
solution.  The  electrical  parameters  are  f = 10  MHz,  a = 10  2 
mho/m.  and  < t ~ 9. 


Tig.  9.  Magnitude  of  the  current  on  the  antenna  shown  in  Kig.  7b. 
I he  interpolation  solution  to  the  problem  is  compared  to  the  exact 
solution.  The  electrical  parameters  are  f = 10  MHz,  a = 10 
mho/m,  and  < F = 9. 
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Kig.  10.  Magnitude  of  the  current  on  the  antenna  shown  in  Kig. 
7c.  The  interpolation  solution  to  the  problem  is  compared  with  the 
exact  solution.  The  electrical  parameters  are  f = 10  MHz,  a = 
10  mho/m,  and  «t  = 9.  Note  that  the  two  solutions  are  in- 
distinguishable. In  order  to  demonstrate  the  effect  of  the  ground  the 
results  from  the  same  antenna  in  free  space  are  also  presented. 
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Kig.  II.  The  Sommerfeld  intergral  fields  must  exactly  cancel  the 
perfect-ground  image  fields  when  the  lower  half-space  has  the  same 
properties  as  the  upper  half-space,  thus  providing  a good  check  on 
the  numerical  results.  Results  are  presented  here  for  the  case  where 
a half-wave  horizontal  dipole  at  a frequency  of  13  MHz  is  lowered 
toward  such  a pseudo-half-space  to  check  the  Sommerfeld  integral 
and  interpolation  portions  of  the  program.  The  input  impedance  is 
seen  to  be  effectively  independent  of  height  to  within  0.003 
wavelengths,  demonstrating  the  consistency  of  the  overall 
numerical  treatment. 
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654.0  ft 


Kig.  12.  Vertical  half-rhombic  used  to  study  the  effect  of  ground  returns.  (The  input  impedance  behavior  over  a perfect  ground  with  and 
without  a ground  return  are  shown  in  the  following  figures,  ground  return  wire  on  interface. ) The  impedance  variations  are  much  greater  for 
the  antenna  having  the  ground-return  wire.  This  is  evidently  due  to  return  wire  and  its  image  acting  as  a transmission  line  to  transform  the  open 
circuit  (with  respect  to  ground)  at  the  termination  end  to  a frequency -dependent  load  at  the  input.  It  might  be  expected  that  over  a lossy 
ground,  this  effect  would  be  much  reduced.  We  have  not  been  able  to  examine  that  possibility  in  this  study,  but  it  is  a topic  that  should  be  pur- 
sued further,  because  of  the  practical  implications  and  to  establish  more  definitively  the  limitations  of  the  numerical  model. 


Frequency  - MHz 


Frequency  — MHz 
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Kig.  13.  Input  resistance  for  antenna  shown  in  Kig.  12.  One  set  of 
data  (a)  Has  a ground-return  wire  and  the  other  set  (h)  does  not. 


Kig.  14.  Input  impedance  for  the  antenna  shown  in  Kig.  12.  on  an 
expanded  frequency  scale.  One  set  of  data  (al  has  a ground-return 
wire  and  the  other  set  (b)  does  not. 


Fig.  15.  Several  sets  of  calculations  were  performed  at  13  MHz 
for  a monopole  antenna  having  various  ground  arrangements.  The 
antenna  and  four  wire  grounds  are  shown  above.  The  input  ini' 
pedance  of  the  antenna  was  studied  as  a function  of  frequency  in 
free  space  and  as  a function  of  height  over  a perfect  ground.  These 
results  are  shown  in  the  following  figures.  Its  resistance  is  seen  to 
vary  smoothly  with  height  between  a free-space  value  of  about  1.55 
12  and  asymptotic  value  near  ground  of  3.20  12.  For  comparison,  the 
resistance  of  the  same  monopole  without  counterpoise  wires  and 
driven  against  ground  is  3.14  12. 


Frequency  — MHz 

Fig.  17.  Input  reactance  of  the  antenna  shown  in  Fig.  15  in  free 
space. 


Frequency  - MHz 

Fig.  16.  Input  resistance  of  the  antenna  shown  in  Fig.  1$  in  free 
space. 


Fig.  18.  Input  resistance  of  the  antenna  shown  in  Fig.  IS  as  a 
function  of  the  height  above  a perfect  ground. 


Fig.  19.  Input  reactance  of  the  antenna  shown  in  Fig.  IS  as  a func- 
tion of  the  height  above  a perfect  ground. 
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\ Fig.  20.  Monopole  antenna  with  (a)  three-wire  and  (b)  two-wire 

counterpoise.  The  inputs  for  the  four-,  three-,  and  two-wire  coun- 
terpoise are  plotted  as  a function  of  height  above  the  perfect  ground 
in  the  following  figures.  The  increase  in  resistance  with  decreasing 
| height  relative  to  the  perfect  ground  case  with  a four-wire  system  is 

indicative  of  increasing  loss  as  the  ground  is  approached.  No  ex- 
planation for  the  bump  at  ~2  ft  in  the  resistance  curve  has  been 
found.  The  resistance  decrease  below  ^0.3  ft  is  apparently  due  to 
numerical  limitations. 

own  radiation  to  that  from  the  antenna  itself,  and  can  be  especially  important  in  influencing  low-elevation- 
angle  coverage. 

Some  results  pertaining  to  the  impedance  effects  of  a radial  ground  screen  and  a vertical  monopole  are  in- 
cluded in  the  previous  section.  Here  we  present  some  examples  of  pattern-control  calculations.  The  antenna  is 
again  a vertical  monopole,  with  tiie  ground-screen  wires.  The  system  was  driven  at  the  points  between  the  an- 
tenna and  screen.  Calculations  were  performed  with  two-wire  ground  systems  as  a function  of  ground-wire 
length  and  angles  between  the  ground  wires.  Note  that  the  presence  of  the  ground  wires  greatly  enhances  the 
radiation  in  a direction  away  from  the  ground  wire  (see  Figs.  23-25). 

This  problem  in  particular  was  explored  to  demonstrate  the  feasibility  of  using  our  approach  to  evaluate 
ground  screen  influence  on  pattern  behavior,  and  not  as  an  extensive  parametric  excerise  per  se.  (It  should  be 
noted  that  a computational  capability  to  do  this  has  not  heretofore  existed.)  While  recognizing  that  more 
aspects  of  pattern  behavior  deserve  attention  (e  g.,  although  we  present  here  only  azimuthal  pattern  results, 
the  elevation  behavior  is  at  least  of  equal  importance),  we  believe  those  results  are  indicative  both  that  the 
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Fig.  21.  Input  reistance  as  a function  of  height  above  a conducting  ground  for  the  antennas  in  Figs.  15  and  20.  The  electrical  parameters  arc  f 
= 13  MHz,  a = 10  mho/m,  and  « r = 9.0.  The  data  for  a monopole  on  a perfectly  conducting  ground  are  also  shown. 
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Fig.  22.  Input  reactance  a,  a function  of  height  above  a con- 
ducting ground  for  the  antenna,  given  in  Fig.  21. 


Fig.  23.  Vertical  monopole  antenna  with  a two-wire  horizontal 
ground  screen  above  a lossy  half-space.  The  two  horizontal  ele- 
ments are  connected  to  the  base  of  the  monopole.  The  electrical 
parameters  are  f = 10  MHz,  j » 10  mho/m,  and  <t  = 4.  The 
radiation  patterns  are  shown  in  Figs.  24  and  25.  Note  that  the  wires 
enhance  the  radiation  very  significantly.  This  increased  radiation  is 
in  the  direction  opposite  to  the  ground  wires. 
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Kg.  24.  Azimuthal  radiation  patterns  at  45°  elevation  angle  for  the  antenna  shown  in  Fig.  23,  when  the  screen-wire  length  is  3 m.  These  data 
are  the  theta  component  of  the  electric  field  and  the  angle  a is  a parameter  as  shown. 


computation  can  indeed  be  performed  and  that  it  provides  invaluable  guidance  for  antenna  design.  For  exam- 
ple, given  a tactical  communications  requirement,  an  approach  like  this  could  optimize  performance  subject 
to  constraints  on  the  numbers,  layout  and  total  ground-wire  length,  antenna  height  and  geometry,  etc. 

Vertical  Half-Rhombic  Parametric  Results 

The  vertical  half-rhombic  (also  known  as  the  lambda  antenna)  includes  as  a special  case  (when  it  is  of 
constant  height)  the  Beverage  antenna.  It  is  essentially  a broadband,  traveling-wave  structure,  of  relatively  low 
efficiency  (5,10%).  Together  with  the  Beverage,  it  exhibits  reasonably  good  forward  directivity  and,  with  ap- 
propriate design  parameters,  can  provide  a low  take-off  angle. 

Calculation  parameters  and  representative  results  are  summarized  in  Figs.  26-28.  (See  also  Table  2.) 
More  complete  results  are  given  in  Appendix  A. 

Sloping  V Parametric  Results 

The  sloping  V antenna  may  actually  be  viewed  as  a modification  of  the  vertical  half-rhombic.  It  is 
derivable  from  the  latter  by  rotating  one  arm  of  the  antenna  towards  the  other  and  feeding  the  resultant  con- 
figuration at  its  apex  rather  than  at  one  ground  connection. 

Calculation  parameters  and  representative  results  for  the  sloping  V are  summarized  in  Figs.  29-35,  with 
more  complete  data  presented  in  Appendix  B.  (See  also  Table  3.) 

A thorough  comparison  of  the  sloping  V and  vertical  half-rhombic  antennas  is  impractical  with  the 
limited  amount  of  data  available  from  these  calculations.  We  can  observe,  however,  that  for  equal  wire  lengths 
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Kig.  25.  Azimuthal  radiation  patterns  at  45°  elevation  angle  for  antennas  in  Kig.  23,  when  the  angle  between  the  two  ground  wires  is  60°. 
These  data  are  the  theta  component  of  the  electric  field  with  the  length  of  the  ground  wires  a parameter. 


and  frequencies  the  vertical  half-rhombic  produces  a maximum  field  strength  roughly  twice  that  of  the  sloping 
V.  The  location  of  their  elevation  pattern  maxima  are  relatively  similar.  Furthermore,  the  vertical  half- 
rhombic  typically  has  a lower  input  resistance. 
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t ig.  26.  V ertical  half-rhombic  antenna  above  a lossy  interface  represented  by  the  x-y  plane.  The  antenna  is  excited  at  a height  H and 
loaded  at  height  H , above  the  interface.  The  total  length  of  the  wire  in  the  antenna  is  I..  A total  of  33  cases  run  for  this  antenna  are  sum- 
marized in  Table  2.  Representative  results  are  shown  on  the  following  figures,  while  the  complete  data  are  given  in  Appendix  A.  The 
tabulated  results  lead  to  three  observations: 

• Increasing  the  length  or  height  of  the  antenna  generally  increases  the  maximum  field  strength. 

• The  maximum  of  the  elevation  pattern  lowers  in  angle  as  frequency  and  antenna  length  increase. 

• The  input  impedances  are  generally  insensitive  to  height  and  decreasing  reactance  with  increasing  frequency. 


Segment  number 


Fig.  27.  Comparison  of  current  magnitude  on  the  vertical  half- 
rhombic  antenna  shown  in  Fig.  26,  when  Hs  = 2.0m  = H,  ,R,  ” 
.100 II,  I.  = 100  m,  f = 10  MHz,  a = 10  ' mho/m,  and  <r  - 4.0  for 
l-VV  input.  The  height  H v Is  marked  by  each  curve. 
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Hg.  28.  Plot  of  the  radiation  pattern  in  the  plane  containing  the  antenna  for  the  same  antenna  and  electrical  parameters  used  in  Fig.  27.  The 
height  H v is  marked  by  each  curve. 
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Table  2.  Results  for  vertical  half-rhombic  shown  in  Figure  26;  ef  = 4,  o = 10  3 mho/m,  RL  = 300  n, 
»s  = Hl  = 2.0  m. 
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fig.  29.  Sloping  V antenna  over  a lossy  interface  represented  by  the  x-y  plane.  The  antenna  is  excited  and  terminated  at  a height  IH  and 
H | above  the  interface,  respectively.  The  elements  of  the  antenna  are  separated  by  an  angle  a.  The  total  length  of  the  antenna  is  I..  A total 
of  22  cases  were  run  for  the  sloping  V antenna,  as  summarized  in  Table  3.  Representative  results  are  shown  on  the  following  figures,  with 
more  complete  data  given  in  Appendix  B.  The  tabulated  data  prompt  three  observations: 

o Increasing  the  wire  length,  the  frequency,  and  the  included  angle  (up  to  %70° ) increases  the  maximum  field  strength, 
o The  maximum  in  the  elevation  pattern  lowers  with  increasing  frequency,  included  angle  and  height, 
o Input  impedance  tends  to  be  primarily  resistive,  except  at  the  lower  frequencies. 
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Current 


t ig.  32.  Radiation  pattern  for  the  antenna  shown  in  tig.  29,  when  f = 10  MHz,  a = HI 
and  R | = 300 II  for  1 -W  input.  The  height  If  f of  each  antenna  is  shown  by  each  curve. 


mho/m.  <t  = 4.  L = 100  m.  H, 
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Table  3.  Results  for  sloping  V antenna,  shown  in  Fig.  29;  = 4.0,  a = 10  3 mho/m,  RL  - 300.0  fi, 

and  Ht  = 2 m. 
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n 

1, 

n 

Voltage, 
at  1 W, 

V 

Power 

in  space 

wave  for 

1-W  input 

0 , 

w 

0, 

W 

Max  field. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  work  described  in  this  report  has  led  to  development  of  a simple,  but  highly  effective  numerical 
method  for  modeling  conducting  objects  located  near  an  interface.  The  procedure  is  based  on  interpolating  a 
“solution  space"  (z  + z',  p space)  between  Sommerfeld-integral  values  calculated  from  straightforward 
numerical  quadrature  or  a series  solution,  to  obtain  the  field  values  required  in  an  integral-equation  formula- 
tion. Because  interpolation  is  a simple  numerical  operation  relative  to  direct  evaluation  of  the  Sommerfeld  in- 
tegrals, a significant  computer-time  savings  can  result,  the  amount  depending  on  the  problem.  For  example, 
using  this  approach  for  modeling  various  types  of  antennas,  we  have  found  a decrease  in  computer  time  of  up 
to  a factor  of  60  compared  with  our  former  method.  The  interpolation  grid  of  course  must  itself  be  set  up  prior 
to  the  integral-equation  solution,  and  this  involves  an  additional  computer-time  investment,  but  once  ob- 
tained, it  can  be  used  for  any  antenna  configuration  that  fits  within  its  boundaries.  Thus,  it  represents  a fixed 
expense  whose  relative  cost  per  antenna  calculation  decreases  in  proportion  to  the  number  of  times  it  is  used. 

An  additional  advantage  of  the  interpolation  procedure  is  the  reduced  probability  that  unnoticed 
numerical  errors  will  invalidate  a calculation;  this  is  because  the  number  of  calls  to  (he  Sommerfeld-integral 
subroutines  is  greatly  decreased.  Furthermore,  after  the  interpolated  field  is  plotted  as  a surface  in  a three- 
dimensional  space,  remaining  errors  can  be  easily  seen,  as  demonstrated  above.  We  emphasize  that,  while  no 
known  errors  remain  in  the  Sommerfeld-integral  portion  of  this  program  (at  least  for  z + z'  > 10  A),  these 

computations  are  too  complex  to  guarantee  that  uniformly  good  accuracy  can  always  be  achieved. 

Another  significant  development  in  this  work  has  been  extending  and  improving  our  numerical 
procedures  to  permit  locating  objects  within  10  '3  X (or  even  less)  of  the  interface.  This  feature  allows  us  to 
model  such  closely  coupled  conductors  as  the  wires  in  a ground  screen,  without  resorting  to  the  approxima- 
tions usually  employed. 
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We  have  demonstrated  the  capabilities  of  the  new  program,  SOM1NT,  by  performing  some  limited 
parametric  calculations  for  several  antennas.  These  initial  applications  are  by  no  means  intended  to  be  con- 
clusive or  complete,  but  rather  to  illustrate  the  kinds  of  information  that  are  derivable  and  the  kinds  of 
problems  that  can  be  handled  with  the  programs.  It  should  be  recognized  that  antennas  located  near  the  earth 
may  experience  perturbations  due  to  a laterally  and  vertically  inhomogeneous  and  undulating  ground  that  the 
idealized  model  employed  here  (a  Hat,  uniform  half-space)  does  not  include.  Whether  such  features  can  ever 
be.  or  even  should  be,  modeled  is  not  obvious.  In  any  case,  the  current  capability  should  allow  antenna 
designers  to  optimize  antenna  performance  while  including  the  dominant  ground  effect  (i.e..  direct  reflection 
from  the  interface)  in  a way  heretofore  neither  analytically  nor  experimentally  possible. 

Most  of  the  effort  in  this  program  has  been  directed  towards  implementing  the  basic  numerical 
procedures.  Much  remains  to  be  done  to  improve  the  program  and  expand  its  applicability.  Some  possibilities 
for  improvement  and  various  applications  are  discussed  below. 


Ground-Screen  Models 

Ground  screens  (or  counterpoises)  are  used  either  to  control  antenna  impedance  or  pattern.  Of  necessity, 
their  analysis  has  been  approximate  and  based  primarily  upon  using  an  effective  surface  impedance.  The  most 
obvious  limitations  of  this  approach  are  its  failure  to  handle  screens  containing  few  or  widely  spaced  wires,  its 
inapplicability  to  predict  low-angle  radiation,  and  its  inability  to  account  for  the  finite  extent  of  the  screen  on 
the  antenna's  current  distribution.  The  surface-impedance  approximation,  furthermore,  does  not  lend  itself  to 
inclusion  in  an  otherwise  rigorous  Sommerfeld-based  approach. 

A really  useful  and  realistic  analysis  of  ground  screens  and  counterpoises  requires,  instead,  that  the 
currents  in  the  ground-screen  wires  be  found  in  a self-consistent  fashion.  This  approach  could  not  have  been 
seriously  considered  without  the  means  to  evaluate  Sommerfeld  integrals  efficiently  and  accurately.  But  the 
techniques  developed  in  the  current  work  make  this  not  only  feasible,  but  practicable.  We  therefore  recom- 
mend that  SOM1NT  be  extended  specifically  to  handle  problems  involving  antennas  that  interact  with  sparse 
systems  of  wires  near  (0.001  A)  the  interface,  and  to  which  they  may  be  attached.  The  use  of  ground  screens 
below  tactical  antennas  could  provide  a significant  improvement  in  their  performance.  With  the  new  method 
described  above,  it  is  now  feasible  to  study  antennas  currently  in  use  to  see  how  ground  screens  might  best  be 
used  or  whether  present  designs  can  be  replaced  by  newer  antennas  with  ground  screens. 


Dielectric  Sheaths 

A recurring  problem,  especially  with  respect  to  buried  conductors,  is  how  a dielectric  or  insulating  sheath 
changes  the  behavior  of  a metallic  object.  This  is  especially  important  in  conducting  media  where  a sheath  pre- 
vents the  direct  flow  of  conduction  current  between  the  object  and  the  medium.  It  is  thus  necessary  to  include 
the  dielectric  coating  in  the  model  if  realistic  results  are  to  be  obtained.  Particular  attention  should  also  be 
paid  to  the  end  conditions  where  the  conduction  current  on  the  wire  may  not  vanish  as  it  does  in  free  space. 
Several  approaches  are  worth  considering  here,  with  perhaps  the  most  promising  one  based  on  a displacement 
current  model  for  the  sheath. 


Two-Half-Space  Problem 

The  current  work  deals  with  wire  objects  located  wholly  on  one  side  or  the  other  of  the  interface,  so  that 
the  Sommerfeld  integrals  are  functions  only  of  z'  (source  height)  + z (observation  height),  and  p (radial 
separation).  Many  applications  occur  where  two  objects  are  on  opposite  sides  of  the  interface.  The  self-fields 
of  these  objects  can  be  handled  efficiently  with  the  present  code.  However,  interaction  of  these  objects  across 
the  interface  involves  Sommerfeld  integrals  of  the  three  space  variables,  z’,  z and  p.  Using  an  interpolation 
scheme  for  this  problem  may  not  be  practicable  because  it  would  then  involve  three  dimensions.  Thus, 
schemes  other  than  interpolation  should  be  considered,  including  the  transmission-coefficient  approximation 
and  various  asymptotic  expansions.  These  and  the  interpolation  technique  could  then  be  compared  to  find  the 
approach  best  suited  to  this  problem  in  terms  of  efficiency  and  accuracy. 
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Effect  of  Nonuniform  Currents 


Associated  with  the  development  of  ground-screen  models  above  is  the  validity  of  the  thin-wire  approx- 
imation as  the  wire  approaches  the  interface.  Specifically,  we  must  assess  the  degree  of  validity  of  the  assump-  * 

tion  that  the  current  is  uniformly  distributed  around  the  wire's  circumference  as  the  wire  approaches  the  inter- 
face. Within  a few  wire  radii  of  the  ground,  this  assumption  appears  questionable,  and  its  effect  on  the  ac- 
curacy of  the  ground-screen  model  is  uncertain.  We  might  expect  that  whether  a thin  wire  is  lying  on  the 
ground  or  a few  radii  above,  the  difference  would  be  small  in  any  case.  That  however,  should  be  studied  to 
validate  the  model.  i 

Penetrating  Wires 

A logical  extension  to  the  two-half-space  problem  above  is  the  case  where  the  object(s)  is  cut  by  the  inter- 
face. This  occurs,  for  example,  when  an  antenna  is  driven  against  a ground  rod  that  serves  as  a counterpoise. 

Many  small  mobile  antennas  used  by  the  Army  currently  employ  ground  stakes  for  connecting  the 
transmitter  to  one  end  of  the  antenna  and  sometimes  for  termination  as  well.  Although  ground  stakes  must  < 

obviously  affect  the  input  impedance  and  radiation  patterns  of  the  antenna,  their  placement  and  design  have  } 

never  been  studied  rigorously.  Computer  modeling  of  this  problem  would  permit  studying  the  effects  of 

ground-staked  antennas  with  the  goals  of  improving  the  present  performance  achievable  with  mobile  I 

antennas.  i 

In  this  case,  we  must  be  concerned  not  only  with  the  field  evaluation  but  also  with  determining  what  j 

special  current  conditions,  if  any,  are  needed  at  the  interface.  For  example,  the  reflection  of  a plane  wave  from  1 

a half-space  requires  continuity  of  tangential  fields  that  can  result  in  discontinuities  in  their  derivatives.  On  a < 

penetrating  conductor,  continuous  tangential  fields  (tangential  H and  normal  E to  the  wire)  at  the  interface  J 

could  result  only  if  the  charge  density  is  discontinuous.  This  suggests  the  need  for  deriving  a current-slope  con- 
dition different  from  that  used  at  junctions  in  free  space.  I 


Antenna  Efficiency 

Previous  work  done  at  LLL  has  demonstrated  that  the  efficiency  of  antennas  like  the  Beverage  and 
horizontal  dipole  can  be  significantly  improved  by  arraying  them.  Efficiency  increases  of  300%  and  more  ap- 
pear readily  achievable,  without  the  use  of  ground  screens.  This  increased  efficiency  has  been  studied  com- 
putationally and  found  to  occur  because  the  multiple  elements  in  an  array  increase  the  rate  at  which  the  an- 
tenna current  radiates,  thus  reducing  the  power  loss  into  the  ground.  Increasing  the  size  of  the  antenna  wires 
has  a similar  effect.  Thus,  antennas  like  the  Beverage,  heretofore  considered  too  inefficient  for  tactical  com- 
munications, might  be  made  competitive  in  this  respect  with  the  vertical  half-rhombic  and  slanted  V,  for  ex- 
ample. This  may  mean  that  ease  of  deployment  and  efficiency  are  not  necessarily  incompatible.  In  this  ap- 
plication of  the  SOM1NT  code,  existing  antennas  as  well  as  alternate  designs  could  be  evaluated  to  determine 
the  extent  to  which  multi-element  arrays  would  be  effective  in  increasing  efficiency,  and  thus  reducing  power 
requirements  of  portable  communications  systems. 


Frequency  Domain  Interpolation 

Although  the  SOMINT  code  operates  nearly  as  fast  as  Antenna  Modeling  Program  (AMP),  which  uses 
the  reflection-coefficient  approximation,  the  cost  of  filling  the  solution-space  interpolation  matrix  is  not  small.  i 

Therefore,  in  evaluating  the  frequency  dependence  of  an  antenna,  it  is  desirable  to  minimize  the  number  of  fre- 
quencies at  which  the  matrix  is  needed.  One  method  is  to  interpolate  between  the  frequencies  at  which  the  | 

matrix  is  calculated.  This  leads,  in  essence,  to  a three-dimensional  interpolation  scheme  whose  storage  require- 
ments could  be  prohibitive.  However,  if  a linear  interpolation  in  frequency  is  used,  only  the  two  solution-space 
interpolation  matrices  that  lie  on  either  side  of  the  desired  frequency  need  be  used.  A solution-space  matrix  at 
the  intermediate  frequency  can  then  be  obtained  from  them;  following  that,  the  procedure  is  identical  to  that 
now  used.  This  new  matrix  could,  of  course,  be  saved  for  subsequent  reuse,  further  reducing  the  required  com- 
puter time.  It  is  worth  noting  that  a similar  approach  could  be  used  for  varying  the  permittivity  and  conduc- 
tivity. 
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Nonwire  Integral  Equations 

The  SOMINT  code  is  currently  implemented  in  an  electric-field,  thin-wide,  integral-equation  code.  But 
many  objects  of  interest  are  surfaces  rather  than  wires.  While  they  may  be  modeled  with  a measure  of  success 
using  wire  grids,  they  can  be  more  efficiently  and  accurately  handled  with  an  integral  equation  applicable  to 
surfaces.  The  same  interpolation  procedure  for  obtaining  the  Sommerfeld  integrals  applies  irrespective  of  the 
form  of  the  integral  equation,  so  that  the  improvements  now  available  in  SOMINT  could  also  be  realized  for 
surfaces. 


User-Oriented  Improvements 

The  Numerical  Electromagnetics  Code  (NEC)  is  probably  the  most  user-oriented  and  thoroughly 
documented  code  available  for  modeling  wire,  surface,  and  hybrid  objects.  It  would  thus  be  very  desirable  to 
add  to  it  the  newly  developed  Sommerfeld-integral  capability  now  included  only  in  SOMINT,  and  to  provide 
the  necessary  documentation.  At  the  same  time,  it  would  be  worthwhile  to  incorporate  further  needed  im- 
provements into  NEC  that  have  been  identified  as  a result  of  our  ongoing  interaction  with  the  AMP-NEC  user 
community.  These  should  include,  but  not  be  limited  to,  improvement  of  the  input  data  generation  routine  to 
minimize  user  errors:  modification  of  the  appropriate  routines  to  eliminate  the  small-argument  errors  that  oc- 
cur when  segment  sizes  less  than  0.001  A are  used;  and  development  of  a routine  to  more  accurately  assess 
impedance-loss  effects  by  integrating  the  loss  term.  Actually,  this  should  be  a modest  continuing  effort  to  re- 
spond to  problems  as  they  arise  and  to  provide  the  information  to  the  user  community. 


User-Experience  Data  Base 

One  problem  that  often  occurs,  especially  when  either  new  users  first  attempt  to  employ  a code  such  as 
NEC,  or  even  when  experienced  NEC  users  apply  it  to  a significantly  different  problem,  is  that  several  itera- 
tions are  required  before  physically  realistic  results  are  obtained.  It  would  be  extremely  beneficial,  by  reducing 
wasted  effort  and  providing  more  timely  results,  to  concisely  document  these  learning  experiences.  If  a 
description;  a given  antenna  type;  the  approach  taken;  results  obtained;  problems  encountered;  etc.,  were  writ- 
ten in  a common  format  and  made  available  to  those  involved  in  antenna  modeling,  greater  efficiency  would 
surely  follow.  As  it  evolves,  this  user-experience  data  base  would  provide  guidance  for  modeling  a growing 
variety  of  applications  and  thus  reduce  the  difficulties  mentioned  above.  It  would  also  provide  guidance  con- 
cerning problem  areas  needing  attention  concerning  computer  code  improvements. 
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APPENDIX  A.  RESULTS  FOR  A VERTICAL 
HALF-RHOMBIC  ANTENNA 


The  vertical  half-rhombic  antenna  is  shown  in  Fig.  A I . All  data  shown  in  this  appendix  are  for  a - 10  "3 
mho/m,  t{  = 4,  RL  = 300  S2,  HT  = 2 m.  Figures  A2-A35  contain  plots  of  the  radiation  pattern  (at  a distance  of 
10 5 meters)  in  a vertical  plane  that  contains  the  antenna  and  the  magnitude  of  the  current  on  the  antenna. 


Fig.  Al.  Vertical  half-rhombic  antenna  abo>e  a lossy  interface  represented  by  the  x-y  plane.  The  antenna  is  excited  by  a 1-W  source  at 
height  H s loaded  at  height  H , , both  abote  the  interface.  The  total  length  of  wire  in  the  antenna  is  L. 
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Current  magnitude 


Hg.  A3.  Results  from  a vertical  half-rhombic  antenna:  (a)  radiation  pattern:  (b)  magnitude  of  the  current  on  the  antenna.  I. 
H,  * 5 m;  freq  = 10  MHt;  R,  OA„  = 300  ii. 
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fig.  A5.  Result*  from  ■ vertical  half-rhombic  antenna:  (a)  radiation  pattern:  (b)  magnitude  of  the  current  on  the  antenna.  L = 100  m; 
H,  = 15  m;  freq  - 10  MHz;  R,  (>Ai>  = J00  Si. 
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Fig.  A6.  Result*  from  a vertical  half-rhombic  antenna:  (a)  radiation  pattern;  (b)  magnitude  of  the  current  on  the  antenna.  I. 
H,  1 2 m;  freq  =*  6 MHz;  R,...,.  *=  300  !2. 
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Fig.  AI2.  Results  from  a vertical  half-rhombic  antenna:  (a)  radiation  pattern:  (b)  magnitude  of  the  current  on  the  antenna.  I.  * 100  m: 
H,  = 5 m;  freq  = 4 MHz;  R,  OAI)  = 300  Si. 
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t ig.  A I t.  Results  from  a vertical  half-rhombic  antenna : ( a ) radiation  pattern ; ( b ) magnitude  of  the  current  on  the  antenna.  I.  = 100  m ; 
H,  = 10  m;  freq  = 4 MHi;  R,  ((A|,  = 300  SI. 
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Fig.  AI4.  Results  from  a vertical  half-rhombic  antenna:  (a)  radiation  pattern:  (b)  magnitude  of  the  current  on  the  antenna.  I.  = 100  m; 
H,  - S m;  freq  - 8 MHz;  Rlm„  - 300  H. 
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Fig.  A 15.  Remits  from  a vertical  half-rhombic  antenna:  (a)  radiation  pattern;  (b)  magnitude  of  the  current  on  the  antenna.  I. 
H.  = 10  m;  freq  = 8 MHz;  R.OAn  = 300  il. 
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Mg.  A20.  Results  from  a vertical  half-rhombic  antenna:  (a)  radiation  pattern;  (b)  magnitude  of  the  current  on  the  antenna.  1.  - SO  m: 
H,  = 2 m;  freq  = 6 MHz;  *|OA|>  = 300 II. 
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Fig.  A2I.  Results  from  a vertical  half-rhombic  antenna:  (a)  radiation  pattern:  (b)  magnitude  of  the  current  on  the  antenna.  I. 
H.  = 5m;freq  = 6MHz;Rln.„  = 300  SI. 


4 


0 0 0?  0 04  0 06  0 08  0.1 


Field  strength  - mV/m 


Fig.  A2  V Results  from  a vertical  half-rhombic  antenaa:  (a)  radiation  pattern;  (h)  mai  nitude  of  the  current  on  the  antenna.  L = 50  m; 
Ht  = 15  m;  freq  * 6 MHz;  R,  ,MI)  * M0  li. 
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fig.  A25.  Rewl  Is  from  I vertical  half-rhombic  antenna:  (a)  radiation  pattern:  (b)  magnitude  of  the  current  on  the  antenna.  I.  = 20(1  m: 
H,  = S m:  freq  = 10  MHz;  R,  OAD  = 300  SI. 
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H,  » 2 m;  fr*q  = 6 MHi;  R1OA0  « 300 !!. 
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fig.  A30.  Results  from  a vertical  half-rhombic  antenna:  (a)  radiation  pattern:  (b)  magnitude  of  the  current  on  the  antenna.  L = 200  m; 
= 10  m;  freq  = 6 MHz;  = 300  S 2. 
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Hg.  \3I.  Results  from  a vertical  hair-rhombic  antenna:  (a)  radiation  pattern:  (b)  magnitude  of  the  current  on  the  antenna.  1.  - 200  m. 
H,  = 15  m;  freq  = 6 MHz:  R(  OAt)  - 300  li. 
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Pig.  A .14.  Results  from  a verticil  half-rhombic  antenna:  (a)  radiation  pattern;  (b)  magnitude  of  the  current  on  the  antenna.  L = 50  m; 
H,  = S m;  freq  = 10  MHi;  R,OA|)  - 300 II. 
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APPENDIX  B.  RESULTS  FOR  A SLOPING  V ANTENNA 


The  sloping  V antenna  is  shown  in  Fig.  Bl.  All  data  shown  in  this  appendix  are  for  a = 10  mho/m,  tr  = 
4,  R(  = 300  il,  and  HL  = Hs  = 2 m.  Figures  B2-B23  contain  plots  of  radiation  patterns  (at  a distance  of  10  s 
meters)  in  a vertical  plane  that  bisects  the  two  halves  of  the  antenna  and  the  magnitude  of  the  current  on  the 
antenna. 


tig.  Bl.  Sloping  V antenna  over  a lossy  interface.  The  interface  is  represented  by  the  x-y  plane.  The  antenna  is  excited  by  a l-W  source 
at  height  H , and  terminated  at  height  H , above  the  interface.  The  elements  of  the  antenna  are  separated  by  an  angle  a.  The  total  length  of 
the  antenna  is  I.. 
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Kig.  B2.  Remits  from  a sloping  V antenna:  (a)  radiation  pattern;  (b)  magnitude  of  the  current  on  (he  antenna.  L = 100m;  H,.  = 10  m; 
freq  = 2 MHz;  a = 60°.  Note  that  here  and  in  the  follow  lag  figures,  the  current  is  shown  on  both  arms  of  the  antenna  with  the  source  in  the 
middle  of  the  diagram. 
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tig.  B4.  Results  from  a sloping  V antenna : (a)  radiation  pattern;  (b)  magnitude  of  the  current  on  the  antenna.  L = 100  m;  H,  = 10  m: 
fre<|  = * MHi;  a « 6fT;  Ra <>ai>  = -tOOU. 
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tig.  B5.  Results  from  a sloping  V antenna:  (a)  radiation  pattern;  ( b ) magnitude  of  the  current  on  the  antenna.  I.  = 100  m;  H(  = 10  m; 
freq  - 8 MHz;  a = 60° : R,  OA„  = 300  Si. 
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Fig.  B6.  Result  from  a sloping  V antenna:  (a)  radiation  pattern:  (bl  magnitude  of  the  current  on  the  antenna.  L = 100  m;  H 
freq  = 10  VI Hi;  a = 60«;  R,(Un  = 300  II. 
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Fig.  B7  Results  from  a sloping  V antenna:  (a)  radiation  pattern;  (b)  magnitude  of  the  current  on  the  antenna.  L = 100  m;  H f = 2 m; 
freq  = 10  MHz;  n - 60° ; R,  OAD  = 300 11. 
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